RIKEN-BNL preprint 



Forming Disoriented Chiral Condensates through Fluctuations 



Dirk H. Rischkc 

RIKEN-BNL Research Center, Physics Department 
Brookhaven National Laboratory, Upton, New York 11973, U.S.A. 



Using the influence functional formalism, classical equations of motion for the O(N) 
model are derived in the presence of a heat bath, in both the symmetric phase as well as the 
phase of spontaneously broken symmetry. The heat bath leads to dissipation and fluctuation 
terms in the classical equations of motion, which are explicitly computed to lowest order in 
perturbation theory. In the broken phase these terms are found to be large for the a field, 
even at zero temperature, due to the decay process a — > nn, while they are small for the tt 
fields at temperatures below T c ~ 160 MeV. It is shown that in large volumes the presence 
of dissipation and fluctuations suppresses the formation of disoriented chiral condensates 
(DCC's). In small volumes, however, fluctuations become sufficiently large to induce the 
formation of DCC's even if chiral symmetry has not been restored in the initial stage of the 
system's evolution. 
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I. INTRODUCTION AND CONCLUSIONS 



At vanishing net-baryon number density and temperatures above T c ~ 160 MeV, lattice calculations of 
quantum chromodynamics (QCD) predict the existence of a phase of nuclear matter where quarks and gluons 
are deconfined and chiral symmetry is restored . One of the primary goals of relativistic heavy- ion physics 
is to create and study this phase in nuclear collisions Q . 

The formation of a so-called disoriented chiral condensate (DCC) has been proposed as a possible signature 
for the restoration of chiral symmetry ||. The idea is the following: in the phase where chiral symmetry 
is restored the quark condensate vanishes, (qq) ~ 0. If at all, in a heavy-ion collision this state can only 
be transiently created. Once the system cools below T c , chiral symmetry is spontaneously broken and the 
system has to evolve back into the true ground state where (qq) ^ 0. If in the course of this evolution the 
pseudoscalar condensate (qTj^q) assumes non- vanishing values (instead of remaining zero, as in the ground 
state), one speaks of a disoriented chiral condensate. 

This mechanism becomes physically most transparent in the framework of the 0(4) model. One identifies 
4>i ~ (qq)i 4>i ~ (9 T i759)j i = 2,3,4, and spontaneously broken symmetry is realized by a potential U(cj>) 
which looks like a (tilted) "mexican hat", with minimum at <f> = (J T ,0) for T = 0. For increasing T, the 
"hat" becomes shallower and the minimum moves towards the origin, such that <fr — > for T — ► T c , and 
chiral symmetry is restored. The creation of DCC's is most likely in the so-called "quench scenario" 0. 
Here, it is assumed that, after restoration of chiral symmetry, the system cools instantaneously to T = 0. 
If the fields are assumed to follow classical equations of motion in the potential U((f>), the evolution of 
the system can be visualized as "rolling down" from the initial state with restored chiral symmetry, (f> = 0, 
into the true ground state cj) — (/ ff ,0). If that happens on a "path" where cj>i ^ 0, i — 2,3,4, the chiral 
condensate becomes "disoriented" . 

The ratio R of neutral pions to the sum of neutral and charged pions was suggested as experimental 
observable p|. If a single domain of DCC is formed, the probability P(R) ~ 1/\R, which is drastically 
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different from the case where pions are emitted in a statistically independent manner, P(R) ~ S(R — 1/3) 
(for large, isospin-symmetric systems). 

The formation and decay of DCC's has been studied in a variety of approaches over the last couple of 
years and the original idea has undergone several refinements. One obvious effect that has an influence 
on DCC formation in heavy-ion collisions is the presence of a background of a multitude of other particles. 
Most of these are pions with typical transverse momenta on the order of a couple of hundred MeV || . In 
it was assumed that these pions constitute a background ("heat bath") of unobserved, thermalized degrees 
of freedom. Their presence leads to temperature-dependent dissipation and fluctuation terms in the classical 
equations of motion which correspond to four-particle interactions between classical fields and particles in 
the heat bath and can be rigorously derived (|] (cf. also |J) via the influence functional formalism j[o|. For a 
typical (average) temperature evolution in a heavy-ion collision it was then studied in 0, how DCC's form in 
an expanding system which is in contact with this (steadily cooling) heat bath of particles. The main result 
was that, on the average, dissipation and fluctuation tend to suppress the formation of DCC's. However, 
fluctuations grow ~ 1/VV for V — > 0. Thus, in small volumes and in a single event, the fluctuations can be 
large enough to destabilize the system and actually enhance the likelihood to form a DCC. 

The authors of Q made two approximations. The first was to compute the dissipation terms in the 
chirally symmetric phase and then to use them for the evolution of the system in the phase where chiral 
symmetry is spontaneously broken. This leads to considerable simplifications, since in the symmetric phase 
these terms are straightforward generalizations of results obtained previously in tfi A theory [S|. For dissipation 
arising from four-particle interactions which are present both in the symmetric and the broken phase, this 
approximation is probably justified in the initial stage of the evolution, where temperatures are high and 
close to T c . It becomes questionable at smaller temperatures due to the fact that not all particles have the 
same mass in the broken phase (the a is heavy, while the 7r's are light). 

There is, however, another reason to reconsider this approximation at temperatures below T c . In the 
broken phase the structure of the underlying Lagrangian is fundamentally different: there are additional 
three-particle interactions. As will be shown in detail in the following, this has the consequence that, 
while the dissipation is ~ A 2 in the symmetric phase and of equal magnitude for all fields c/) a , a = 1, . . . , 4, 
dissipative corrections arise already to first order in A in the broken phase, and are of sizable magnitude 
for the a degree of freedom and rather small for pions. The former correspond physically to the decay of a 
a into two 7r's (cf. also jy]), and are non-vanishing even at T = 0. Therefore, the dissipation coefficients 
below T c are different from an extrapolation of the results obtained in the symmetric phase. 

The second approximation made in J?j] was to infer the variance of the fluctuation terms £ a , a — 1, . . . , 4, 
from the dissipation coefficient rj via 

(Ut)Zb(t')) = ^5(t-t')S ab . (1) 

Apart from the fact that this equation does not account for different dissipation coefficients for a and 7r's, 
one has to note that the factor 2T stems from the high-temperature (i.e. classical) limit of a more general 
expression (cf. M and below). For temperatures T smaller than the typical mass scale of the theory, i.e., for 
T < m v ~ T c , one therefore expects sizable deviations from (^). In particular, whenever rj happens to be 
finite at T = (which, as mentioned above and shown below, is indeed the case for the a field in the broken 
phase), eq. ([!]) predicts that fluctuations vanish even in the presence of dissipation, in contradiction to the 
dissipation-fluctuation theorem. 

The aim of this paper is to make a first step towards a consistent treatment of dissipation and fluctuation 
in the framework of the 0(4) model in the phase where chiral symmetry is broken. The outline of the paper, 
as well as the main results and conclusions are as follows. In section II, the derivation of the influence 
functional is presented for a system of N real-valued scalar fields <f>\, . . . , 4>n- Short wavelength modes, i.e., 
those with "hard" momenta |k| > k Cl where k c is an arbitrary momentum scale, are separated from long 
wavelength modes, i.e., those with "soft" momenta |k| < k c [Q, and the reduced density matrix for the soft 
modes is obtained from the full density matrix by tracing over the hard degrees of freedom, which thus are 
thought to constitute the unobserved background ("heat bath") mentioned above. The influence functional 
enters as a phase factor in the reduced density matrix. This section is a straightforward generalization of 
the treatment in H . It serves merely to introduce the notation, and can be skipped by readers familiar with 
the subject. 

Section III contains the derivation of the classical equations of motion for the soft fields by expanding 
the reduced density matrix around its diagonal elements. The main difference as compared to previous 
treatments of the subject |q,0,n3| is that functional derivatives of the influence functional are expressed as 
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averages over functional derivatives of the action characterizing the interaction between soft, classical fields 
and the unobserved hard degrees of freedom. This has the advantage that, in a perturbative computation of 
the correction terms to the classical equations of motion up to some given order n in the coupling constant, 
expectation values of the hard degrees of freedom have to computed only to order n — 1 . It is furthermore 
shown that for systems with more than one field degree of freedom, there can be cross correlations between 
the noise terms for different fields. 

In Section IV, the general framework derived in the previous sections is applied to the O(N) model in the 
symmetric case. The corrections to the classical equation of motion are first computed to order A, where 
they only change the mass term in the classical equation of motion, and then also to order A 2 , where they 
lead to dissipation and fluctuations. The treatment is fairly cursory, since this case is rather similar to (j) 4 
theory (or, in other words, the 0(1) case) discussed in detail in ||. The main focus is to demonstrate the 
applicability of the method developed in Section III. It is shown that the damping coefficient agrees with 
previous results in the cases N = l H and TV = 4 Q. 

In Section V the O(N) model is discussed for the case of spontaneously broken symmetry. In this case, 
only corrections up to first order in the coupling constant are considered, but due to the presence of two 
interaction vertices in the Lagrangian (one proportional to A and the other proportional to A/^ ~ A 1 / 2 ), 
there is dissipation and fluctuation already to this order in A. In particular, the dissipation coefficient for 
the a field is shown to be large even at T — 0, r\ a ~ ra a ~ 600 MeV. The physical process responsible for 
this is the decay of a a into two 7r's. This has important consequences which are discussed in Section VI. 
On the other hand, the dissipation coefficient for the 7r fields is small for the temperature range of interest, 
the reason being that scattering of a ix or u from the heat bath off a classical it field is strongly suppressed 
by phase space. For T — > 0, as well as in the chiral limit — > for arbitrary T, one even has tj^ — > 0. The 
classical equations of motion for er's and 7r's are derived and studied in detail for the k = modes of the 
classical fields. It is shown that the static solution for the a field corresponds to the well-known shift of the 
vacuum ground state at finite temperature. Moreover, the validity of Goldstone's theorem at the classical 
level is checked. 

In Section VI, arguments are presented that the formation of DCC's is most likely in a quench scenario, 
i.e., at T = 0. Numerical solutions of the equations of motion for the homogeneous modes of the fields at 
zero temperature are then presented, which show that, in large volumes, the large dissipation coefficient 
for the a field leads to a rapid damping of oscillations of all classical fields, including the pions. The 
formation of DCC's seems thus not very likely in large systems. In small volumes, on the other hand, the 
fluctuations associated with the dissipation are large enough to disorient the pion fields and possibly lead to 
the formation of DCC's. These results are in agreement with those found in Q], except that here they apply 
even at T = 0. This may have the experimentally interesting implication that DCC's are perhaps formed 
more readily in collisions of lighter ions, or even in pp-collisions (which constitutes a possible explanation 
for the CENTAURO events |15|), while they are presumably less likely to be formed in collisions of heavy 
ions. A definite conclusion, however, can only be drawn after performing dynamical simulations including 
modes with finite k and taking the overall expansion of the system into account fl6| . 

There is, however, another possible consequence of the results found here. As long as the volume is 
small, ~ 10 fm 3 , the disorientation of the 7r fields (and possibly DCC formation) becomes likely even if the 
system's evolution starts near the true ground state (even small perturbations in the 7r fields suffice). As a 
consequence, in small systems restoration of chiral symmetry does not seem to be a necessary prerequisite 
to observe disorientation of classical pion fields. It is, however, rather likely that this fluctuation-induced 
phenomenon is related to ordinary fluctuations in finite volumes. 

The results of the present work have to be viewed in the light of the following two comments: (a) the 
classical approximation works well in the limit of large occupation numbers. For instance, in thermody- 
namical equilibrium this is achieved for modes with energy ui <C T. However, the lowest-energy mode for 
(non-interacting) a particles has to — m„ ~ 600 MeV which is much larger than the temperature in the 
broken phase (T < T c ~ 160 MeV). Therefore, at least in thermodynamical equilibrium, the a field should 
not seriously be considered classically. In a sense, the expected large quantum corrections become manifest 
in the fluctuations induced by the decay a — > tttt as calculated in the present work, (b) The dissipation and 
fluctuation terms are here computed to first order in perturbation theory. However, for realistic parameters 
of the O(N) model the coupling constant A ~ 20, which renders the perturbative expansion uncontrollable. 
Future studies will have to improve on this point. 

Possible other extensions of the present work are (a) the inclusion of collisional interactions between cr's 
and 7r's which are of order ~ A 2 |17|| , (b) the study of long wavelength modes with finite momenta k instead 
of the homogeneous modes only, as well as (c) the study of DCC formation including these effects with a 
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realistic temperature evolution |lq| . 

Units are % = c = fee = 1, and the metric tensor is g^ v = diag(+, — , — , — ). 

II. DERIVATION OF THE INFLUENCE FUNCTIONAL 

Let us consider a quantum system of N real-valued scalar fields (j> — (<f>x, . . . , </>at), characterised by a 

Lagrangian density C{<p). The Hamilton operator H = J d 3 x H(<j>(x)) , Ti = w ■ d t (f> — £, tt_ = dC/d(d t (j)), is 
assumed to have no explicit time dependence. The time evolution of the density matrix reads: 

id t p=[H,p], (2) 

with initial condition p(U) = pi. The formal solution at time tf is: 

Kt f ) = U(t f ,t i )p i U(t i ,tf) , (3) 

where 

U(tf,ti) = exp{-iH(tf-ti)} = U\U,tf) = U- X {t u t f ) (4) 

is the time evolution operator. Let us choose a basis |l<^)| °f eigenfunctions of the Schrodinger field 
operator ^.(x). In this basis, one finds for the density matrix element 

p{± r 4L f \t f ) = ^ f \p{tf)\^ f ) = J D4.D£p(£.,£;t i ) J~ f V$J~ f V$ exp{i(S[@- S[£})} , (5) 

where one has employed eq. (|^), the completeness relation 

1 = J D£ |<£.) <£.| (6) 
(here D(f) = J\a=i IIx d0a(x)), and the path integral representation 

(0 | U(t,,U) |£) = VI exp {i S\$} , (7) 

where = ]\ * n t)X d<p a (t, x) , and S\£ = j t [ f dt J d 3 x£(0(i,x)) = J t [ f & A xC{^{x)). Note that it is 
customary ||,[l3|,[l4| to employ the closed-time-path formalism ([LSI to simplify the right-hand side of eq. (||) . 
Although this is an elegant bookkeeping device, for the sake of clarity we shall continue to work in terms of 
the fields 4> and </>'. 

At each time t, a particular field component with 3-momentum k has the Fourier representation 

0(f,k) = / d 3 xe- ltx ^(t,x) . (8) 



Let us now separate "hard" from "soft" degrees of freedom [p|Jl2j| . More precisely, let us define soft fields 

/d 3 k 
-_e ik -£(t,k)e(fc c -|k|) , (9) 

and hard fields 

£(t,x) = J ^e lk x 0(t,k)9(|k|-fc c ) , (10) 

where k c is an arbitrary momentum scale separating hard from soft momentum modes. Obviously, the space 
spanned by {|^>)} is the product space spanned by {\f}} and {|$}}, {]<!>)} = {If)} ® {!$)}■ (This is most 
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easily seen in the space of functions (j>(t, k).) Also, 0(x) = p>(x) + and 5[0] — S[p] + + Si[(p,^\, 

where Si is the action characterizing interactions between soft and hard fields. 
The density matrix (^) then assumes the form 

p(<P f ,*f,<F! r &f',tf) = J D£ D$, D£ D*J p(^ , ; ^ , ) 

x f"' P<p f~ f P$ 2V /" _/ exp {t (,%] + S[$] + Si[tp, t] - .%'] - 5[*'] - 5/[^', $'])} . (11) 



Let us now assume that the interactions between soft and hard fields vanish at the initial time ti. Then, the 
initial density matrix is block-diagonal, p, = pf^ eg) and 

p{&&\£&\k) = P {v) (& ) &U)p { * ) &,$bU) . (12) 

The object of interest in the following is the reduced density matrix p^{tf , tp' ;tf) for the soft fields, which 

is obtained by tracing (|ll| ) over the degrees of freedom of the hard fields. This reduced density matrix has 
the form 



where the influence functional Sif[<P,<p'] is defined by JlO| 

exp {iSi F [ip,(/]} 

r f± s r± t 

= / D^DtjDtJ/jWdj.Sj;**) / X>* / X>*' exp{^( > S'[t] + S'/[<£,$]-S'[$ , ]-5* 7 [^ , ,$ , ])} . (14) 

J 7$. 

An obvious property of the influence functional (for real- valued scalar fields) is 

SifYPj^] = ~S* IF WM • (15) 
For the proof, note that for any element of the (hermitean) density matrix one has p*(0, 0':t) = 

(0| p(t) = m P f (t) 10) = m m 10) = pW, 0; *)■ 

Another property is 

Si F [vM=U. (16) 

For the proof, note that in Sif, p plays the role of an external, and thus fixed, field. Defining the action 
S^j^?] = + Sj-[</j,$] for hard fields in the presence of this external field ip, and a corresponding time 
evolution operator U v (tf,ti) with matrix elements 

{* f \U v (tf,ti)\*i) = /"~'z>£exp{iS p [$]} , (17) 
J* 

one then reverts in eq. ( |l4| ) the steps which led to eq. (|5|). Then, exp {z Sif[<P, ¥>]} = / p^f' ($^, tf) = 
Tip^p\tf) = 1, consequently, Sjir[<p,<p] = 0. Here, p^p (tf) is the density matrix for the hard degrees of 
freedom as it evolved from its initial value pffi to the final time tf subject to the time evolution operator 
Upitfjti). This evolution in general differs from the one when the external field ip is absent. 



III. DERIVATION OF CLASSICAL EQUATIONS OF MOTION FOR THE SOFT FIELDS 

The classical equations of motion for the soft fields <p are determined by expanding the phase S[tp] — S[ y/\ + 
Sif [<p, p'} in eq. ( |l3| ) around field configurations (p = cp', i.e., where this phase vanishes according to (|l6|) . 
Let us introduce 
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£ = 9 (£ + £') ' A £ = £ ~ £' ' 



(18) 



and expand 5[</j] — 5[</s'] 4- Sif[^P, f'] to quadratic order in A<y9: 

5M - 5k/] + 5/Hrf 

J5kp] 1 fSSiF [<p,(p'] 6Si F [<p,<p'] 



d^x 



Stp(x) 2 \ Sip(x) 



5<p'(x) 



Aip(x) 



i r*f 



d 4 xd 4 y 



(19) 



PSi F [<p,<f?\ PSi F [<p,<pf\ PSi F [<p,<ff] PSi F [<p,<ff\ 



a, 6=1 



8<p a (x) 6<p b (y) 8tp a (x) 5ip' b (y) Stp' a (x) 5p b (y) Sip' a (x) 5ip' b (y) 



<P=<P =v 



Previous derivations of classical equations of motion in a background of hard fields usually compute the 
functional derivatives of the influence functional in eq. ( |l9| ) directly for the specific system under consider- 
ation. However, these functional derivatives can be further evaluated also in the general case. To this end, 
let us define an "n-point function": 

Aato^^K^i)^^)---^^^)^ 1 ^)^ 2 ^)---^^) > » = > (2°) 

i=l J=l 

of the fields $, Here, aj and 6j label components of these TV-dimensional fields. Without loss of generality 
one may assume that t\ > t% > ■■■ and t' x > t' 2 > ■■■■ Let us then define the average of A ($,$') in the 
presence of the "background field " (p as: 

x / D$ / D$' exp {i (5[t] + 5/[<£, $] - 5[$'] - 5j[g, $'])} A($, $') . (21) 

Using the time evolution operator ( p"^ ) (with ip = cp), one notices that this is equivalent to 

(A($,g))e = Tr{$\t f )A&g)} , (22) 

i.e., the usual expectation value of the n-point function A in the (in general, non-equilibrium) ensemble 
characterized by the density matrix p^\tf). Note that in this expectation value all <& fields are time- 
ordered, while all J>' fields are anti-time- ordered. 

Writing Sif = —i In exp{i Sif}, one now observes with the definition (14) and the property ( |l6| ) of the 
influence functional, and the definition of the average @ that 



5Si F [<p,<p'] 



Stp(x) 
SSif [<p,<p'\ 



5tp'(x) 
5 2 SifVp, ¥>'] 



6ip a (x)6(p b (y) 



<p=<p'=<p 



ip — if' — if 



S(p a (x)S(p' b (y) 



SSi[£,$\ \ 
<*£0) /_ ' 

/SSi\^g]\ 

\ 5 ^ x) l 9 ' 

' 5 2 S 1 [^] \ 
5(p a (x)8(p b (y) / 

5<p a (x) 5(p h {y) 
1 5Si[(p,$\ 5Si[tp,& 



(23a) 
(23b) 



SSi[(p,m /tf5i[^,*] 



S( Pb(y) 



5(p a (x) 5(p b (y) 



+ i 



(23c) 
(23d) 



G 



5 2 S IF [<p, l p'} 



S(p' a (x)S(p' b (y) 



_ \Sipa(x)Sif b (y) J _ 

p — p'—p \ I p 

J SSifeg] SSi\^g] \ _./SSi\^g]\ /6Si\&g]\ 

\ Sif a (x) 6<p b (y) I \ Sip a {x) I \ 5ip b (y) ' 



Since the fields <p, <E>, are assumed to be real-valued, eq. ( |l9| ) can be further simplified using symmetry 
properties of the average with respect to exchange of J> <-> <j?', for instance: 

/ 6Si\£,§L'] \ _/ SS I [m \ (24) 



\ Sip(x) / \ 8ip(x, , 



Then, eq. ( pj| ) becomes: 

d 4 xd A y 

,.b=l 



ft} i rtf JL 

S[<£\ - S\g\ + S IF \<p,</\ ^ / d 4 x£{x) ■ A<p(x) + - / d 4 xd 4 y V A<p a (x)l ab (x,y) A<p b (y) , (25) 



where 



8S[(p) 

£(x) = T= f{+n(x), (26a) 

0(p(X) 

K{x) = Re ( ;_7 ~ J ) , (26b) 



Stp(x) 

Z a b{x, y) = ~ Im ( F= ? ^ ? \ ) - Re ( - > Re - 



2 \ 6<p a (x) 6<p b (y) I \ 5tp a {x) / ^ \ 5ip b (y) , 9 



4 \ \ &?a(x) S(p a (x) J y 8<fb(y) 5(p b {y) 



(26c) 

are real-valued functions. Inserting this into expression ( |l3|) and changing the integration variables tp, tp' to 
tp, Atp, one obtains 



A £/ ; i/) ~ y Dg. DA^. (g., A^; ti) y_ 2?£ VA^ 

xexpji^" d 4 x£ (x) ■ Atp(x) - i jf d 4 £ d 4 y A<p a (x) l a b{x, y) Aip b (y) } . (27) 
The term quadratic in Ay> in the argument of the exponential can be rewritten introducing auxiliary fields 

i=(a,...,6v), 

exp f d 4 xd% Atp a (x)l a b(x,y)Atp b (y)\ = J 2?£P[£,£] exp ji jf d 4 a;£(a:) • A^(x)| , (28) 
where 

P[£,g] = A/"|£] exp j-5 d 4 xd 4 y £ kto^fotfj&fo) j ( 29 ) 

is a normalized Gaussian measure. Note that P in general depends on tp through I. It is now possible to 
perform the functional integration over Aip (except for the one at t{) in expression (B^), with the result 
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pM(£ f) A£ f ;t f ) ~ JviJ Dg. DAp. pM(g., A £ .; U) J f Vp P\£,@ S [£(x) + £(»)] . (30) 

This result means that, for times ti < t < t/, the functional 5 function forces the fields Cp to obey the classical 
equations of motion 

5S\0\ 

-£( x ) = -^L-1l{x)=Z(x) . (31) 

0Lp(X) - 

The condition —5S[0\/5ip(x) — is the usual classical equation of motion. The new term TZ characterizes 
the interactions of the soft, classical fields <p with the hard, unobserved degrees of freedom. As it will 
become clear in the following, a part of these interactions is to be interpreted as dissipation. The associated 



fluctuating noise field is represented by £ on the right-hand side of the equations of motion (31). These 
equations of motion are therefore Langevin-type equations. Note that in general the noise is not white, 
since X need not be proportional to 8(xq — yo), and can be multiplicative, due to the dependence of X on (p 
HfL4[ . It should be mentioned at this point that multiplicative noise terms are treated slightly differently in 
than in the present work. While here there is only one noise field with a (^-dependent variance, 
different noise fields are introduced with ^-independent variances. 
Another interesting aspect for systems with more than one field degree of freedom is that, since X is 
in general not diagonal, X ab (x,y) ^ for a =/= b, there can be correlations between the noise terms in the 
equations of motion for two different field components (p a and (p b . In the following, this general formalism will 
be applied to the O(N) model, both in the symmetric case as well as with spontaneously broken symmetry. 




IV. THE O(N) MODEL IN THE SYMMETRIC CASE 



The Lagrangian of the O(N) model is: 



1 



C{fa) = ^{d ll fad^-rrv l fa 



A_ 

N 



(32) 



where <f> = (fa, fa, ... , 4>n) is an ./V-dimensional vector of real- valued scalar fields. 

In the symmetric case, m 2 > 0. The vacuum state of the model is at the minimum of the potential 



U(fa 



A 

TV 



(33) 



which is vac = 0. Decomposing the field <f> according to eqs. (||,[lo|), the action corresponding to interactions 
between soft and hard fields reads: 



2A 

' AT 



d 4 3 



2ip(x) ■ (p(x) <£_{x) ■ $(x) + (p(x) ■ ip(x) $(x) • + 2 \cp(x) ■ M.x)\ 

+ 2tp(x) • • . (34) 



Consequently, 



4A 
~N 



<5 2 ff/[^,$] 

5<p a (x)S<p b (y) 



4A 



+ 2£(a0 • ($(x) $ a (x)) 9 + • t(x) * a (x)) 9 



(2 6 ab ftx) ■ ($(x)) £ + S ab (Mx) ■ £(a;)> £ + 2 (p a {x) 
+ 2 (p b (x) (**{x)) 9 + 2 ($ Q (x) * b {x)) 9 ) 6^(x - y) . 



(35a) 



(35b) 



In the following, the averages on the right-hand side will be computed perturbatively in the coupling constant 
A. The expansion in powers of A of the exponential in the integrand of eq. (|2l|) reads: 
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exp {i (Sm + Si[v,$\ - - Szil,®'})} ~ exp {i (S„[£] - [£'])} 

x l-^| f d 4 j/ (4£(y) • g(2/) £(») • [£(») - f(t/)] + 2£fo) • £(») [£(y) • $(y) - £'(y) • f(y)] 

+ 4 [£(y) • £(y)] 2 - 4 [g(y) • £'(y)] 2 + 4£(y) • [$(») *(») ■ l(y) - £'(y) £'(») ■ !'(</)] 
+ ■ $(y)] 2 - [*'(») ■ ] 2 ) + 0(A 2 )] , (36) 

where So[$] is the action for non-interacting hard fields. The quartic self-interaction term of the hard fields 
has been included in the perturbative treatment of Si. Let us define an average ( • ) in analogy to (|2l|), 
where S[$] + Si[<p, $] - S\&] - S/ [£,$'] is replaced by So[<5] - So[$']- If we additionally assume the initial 

density matrix p^ to be of the form p^ = exp{— Hq/T}/Z, this average is then the usual thermal average 
in a non-interacting system at temperature T. Since this average involves a Gaussian measure in function 
space, the average (|2l| ) of an arbitrary n-point function -A ($,$') vanishes for odd n and can be decomposed 
into a sum over products of 2-point functions for even n. These 2-point functions are p9[ : 

(*„(a;) Mv))o = Tr {pW T ($ a (x) My)) } = S ab D++(x - y) 

= 8ab[D > {x-y)Q{x -y ) + D < {x-y)Q{y ~xo)] , (37a) 
($„(a!) *i(»))o s Tr <l£(y)} = <U C < (x - y) , (37b) 

«(x)$ fc (y)) = Tr {^^(^(y)} =*„6 D>(a-i/) , (37c) 

(*' a {x) %{y)) = Tr {pW f ($;(x) $' 6 (y)) } = <U D-(x - y) 

= S ab [D < (x ~ y)Q(x ~ y ) + D y (x ~ y)e(y - x )} , (37d) 

where T stands for time ordering, T for anti-time ordering, and where translational invariance in space-time 
has been assumed. The functions D > and D < have the Fourier representation 

A(*,x) = J ^9(|k|~fc c ) e ikx A(i,k) , i = > or < , (38) 

where 

D > (t,k) = JL{[l + n (£ k )] e -'^' +Ii (£ k )e ,£kt } , £><(*, k) = £>>(-*, k) . (39) 
Here, u(x) = (e x ^ T — 1) _1 is the Bose-Einstein distribution function, and = (k 2 + m 2 ) 1 / 2 . 



A. Interactions between soft and hard fields to first order in A 

To determine the interactio n ter m TZ, eq. ( [26b| ), to order A, one needs to compute the expectation values 
on the right-hand side of eq. ( 35a ) only to order A = 1 , since there is already an overall factor of A on the 
right-hand side. This means that the averages ( • )» on the right-hand side can be replaced by ( ■ )o. This 
simplification can be traced back to the fact that the functional derivatives of the influence functional in eq. 
( |l9| ) were rewritten in terms of averages over functional derivatives of Si, cf. eq. (p3|), which itself is already 
of order A. Evaluation of the averages therefore yields 

<$a(a%^ ($a(s)>o =0 , (40a) 
($„(x) $ 6 (s)>£ -» ($«(») $i(aO>o = S ab D ++ {0) , (40b) 
($ a (i) $ 2 (x)) £ -» ($ a (x) $ 2 (x)) - . (40c) 



Thus, 



n a {x) = - 4 {N + 2) A Re D ++ (0) y a (z) . (41) 
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With eqs. (37a), 



and (|3^) one derives: 
^++(0) ' 



e(|k|-fc c ) 



( 27r )3-v-- -- 2 £ k 



[l + 2n(£ k )] 



(42) 



This expression diverges, due to the vacuum contribution corresponding to the 1 in brackets. The diver- 
gence can be removed in the standard way, for instance by introducing an appropriate counterterm in the 
Lagrangian. In the following, it is implied that, wherever necessary, such divergences have been removed 
accordingly. Note that for massless particles and in the limit k c — ► 0: 



d 3 k 1 



(2tt)3 E k 



11 
12 



(43) 



Since the expression ( |42| ) is real- valued, one obtains for the left-hand side of the classical equation of motion 
for the field component Cp a : 



£ a (x) 



□ +m 2 (T) + i^^(x)-(£(x) 



Here. 



m 2 (T) 



4 (iV + 2) A 
N 



£++(0) 



(44) 



(45) 



is the usual thermal mass (squared) to first order in A. The term \m 2 {T) — m 2 ](p a (x) can be graphically 
depicted as in Fig. |]. In this and the following graphs, external legs correspond to classical fields (p a , while 
internal lines correspond to propagators -D++ of hard modes. The thin external leg on the left side of 
the vertex has no correspondence in the classical equation of motion, it corresponds to a factor Aif a (x) 
multiplying E a {x) in the argument of the exponential function in eq. (g7|). It is included here to show that 
the interaction in principle involves four particles. 




FIG. 1. The order A contribution to the thermal mass. 



The result ( f44| ) means that to first order in A, the only effect of the presence of the hard modes (the heat 
bath) is a modification of the mass term in the classical equation of motion for the soft fields. There is 
obviously no dissipation to this order in A. Consequently, the fluctuating field £ a has to va nish , too. This is 
verified by explicitly computing X a b(x,y). One first notes that the last two terms in eq. (26c) are of order 
A 2 and thus can be neglected to first order in A. With eq. (35b), the remaining term yields: 



1ab{,x,y) = S a b 



4 (TV + 2) A 
2N 



ImD ++ (0)(5 (4) (a;-y) 



(46) 



which obviously vanishes, since the expression (fl2| ) is real-valued. Therefore, to first order in A, the intro- 
duction of a fluctuating noise field £ via ( p8| ) is obsolete. This also implies that, to this lowest order in A, 
the classical field (p{x) has no means to equilibrate with the background of hard modes. One has to go to the 
next higher order in A to achieve this. 



10 



B. Interactions between soft and hard fields to second order in A 



To second order in A, the expectation values on the right-hand side of (|3^) have to be evaluated to order 
A. The result (neglecting terms which vanish on account of momentum conservation) is: 



(*«*(«))£ 



-i-jy I d 4 y [D {■>■- i/)-D U ■ //)] £l//) ■ £!//) -„(//) . 



S ab D ++ (Q) - i ^ J ° d 4 y [Dl(x -y)-D%(x- y)\ 

x S ab £(y) • £(y) + 2 ^ a (y) £ 6 (y) + <U (N + 2) £>++(0) 
$ (z)<I> 2 (x)) £ = -i^(l + 2<U)^ °d 4 y [D»(x-y)-D»(x-y)] ^(y) . 



(47a) 

(47b) 
(47c) 



In eq. (47c) a term was omitted which vanishes on account of momentum conservation in the expression 
/ d 4 x£_ (x) ■ Aip(x). The left-hand side of the equation of motion (EHI) is then: 



£a{x) 



a + m 2 {T) + ^l(x)-l(x) 



2 3 



where 

^(x) 

T a (3) (x) 



d 4 y [X>>(x - y) - D < (x - y)] \tp(x) ■ <p(x) cp a (y) + 2(p a {x) tp(x) ■ tp(y)] <p(y) ■ <p(y) 

/wo 

d 4 y[Dl(x-y)-D 2 < (x-y)] [(N + 4) <p a (x) £(») • £(y) + A(p(x) • g(t/) £ tt (y)] , 

/SCO 

d 4 j/ [£>|(x - y) - D*(z - j,)] 2(iV + 2)^(y) , 



and where the thermal mass m(T) is given by: 

.4 (TV + 2) A 



2 (T , 2 , 4(iV+2)A 
m (1 ) — m H — 



jH d 4 y [D 2 (ar - y) - £>* (x - y)) ) £ ++ 



(0) 



(48) 

(49a) 
(49b) 
(49c) 

(50) 



Ineqs. ([48|) - (|50|) , use has been made of the fact that D% [x—y] = D > (x— y), and thus -D™ (x— y)— D™ [x— y) = 
2ilm_D™ (x — y), such that all expressions in these equations are real-valued. For N = 1 and A — > ,g 2 /4!, eqs. 

( |47| ) - ( |50| ) reduce to the corresponding ones found in || in tfi 4 theory. The terms T a and the order A 2 thermal 
mass correction are graphically depicted in Fig. |2| Since — i[D™(x — y) — D<(x — y)] = 2ImDV > (x — y) = 
2 Im DV: , (x — y) for xq > yo, in order to make this graphical correspondence it is implied that one has to take 
the imaginary part of any combination of propagators linking the space-time points x and y in the diagrams 
of Fig. |2| (as well as in the following figures). Since all external legs correspond to real- valued classical fields, 
however, this is equivalent to taking the imaginary part of the whole diagram. Note here the advantage of 
working with the functions D instead of G = i D ]q,|l4|| , which would require taking both real and imaginary 
parts of the diagrams, depending on the number of internal lines. 

The two graphs (a,b) correspond to the two terms in the expression (49a). A su m ove r the external legs 
b and c is implied. Graphs (c,d) correspond to the two terms in the integrand in (49b). A sum over the 
external legs b is implied. Figure || (e) corresponds to expression (49c). Finally, Fig. £ (f) is the second order 
contribution to the thermal mass (times (p a (x)). 
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(e) 




(f) 



FIG. 2. (a,b) The two terms in . A sum over b, c is implied. (c,d) The two terms in 7^ 2 ' . A sum over b is 
implied, (e) The term Ta 3 \ (f) The order A 2 contribution to the thermal mass. 

is Fourier 



In order to study the time evolution of the classical momentum mode functions <p a (t, k), eq. (|4 
transformed with respect to x to obtain the classical equation of motion for these modes: 



4 A f d 3 pd 3 q 



4 A 



2 3 



+ 



N 



(2tt) 6 



Q(k c - |p|)6(fe c - |q|)6(fe c - |k - p - q|) (p{t, p) • <p(t, q) <p a {t, k - p - q) = £ B (t, k). (51) 



Here = [k 2 + m 2 (T)] 1 / 2 . In order to see how dissipation enters in the equation of motion ([sil), one has 

to compute the three terms Ta l \t, k), i — 1,2,3. This calculation is rather similar to the one presented in 
H for 4 theory. The details are therefore referred to Appendix A and only the main results are outlined 
here. From eq. ([19]) it is obvious that the time evolution of the classical fields prior to time t influences the 
value of the field at time t. The equation of motion (|5l]) is therefore a so-called delay -differential equation. 
While such integro-differential equations can in principle be solved, it is much simpler to make an Ansatz 
for the time evolution of the fields (p(y) in the terms (49), which renders the equation of motion local in 
time. Such an Ansatz is the so-called linear harmonic approximation [g], cf. Appendix A. In particular, for 

the term Ta {t, k) one then obtains 

T a ^(t,k) = -2 (Ar + 2)6(fe c -|k|) 



du M3(uj, k) 

27T Ek — LO 



<fa{t,\a) 



■M 3 (£k,k) 
2£ k 



d t tp a (t, k) 



where P stands for the principal value and A^w, k) is the Fourier transform of A4s(x) = D^(a 



(52) 



DUx), 



cf. eq. (A2c). The first term represents a thermal mass correction for the field <ys a (£,k), while the second 
yields a damping term +rj(k) dt(p a {t, k) in the equation of motion (|5l]), where the damping coefficient is: 



tiQl) 



4 A 
~N 



2{N + 2) 



■M 3 (£k,k) 
2£ k 



(53) 



For N =1 and A - 
(cf. eq. (41) of |). 



g 2 /4!, this agrees with the result (52) of ||], noting that their function iAi^ = g^M^/Q 
For k = and k c — * 0, the evaluation of .M3 simplifies (cf. Appendix B) and yields 
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where Li2(x) is the dilogarithm (or Spence's integral, cf. eq. (JB7|)) . For (j) 4 theory, i.e., N = 1 and A — ► <? 2 /4! 
the damping coefficient is 

which is twice the on-shell damping rate at zero momentum, computed to 2-loop order in </> 4 theory |2(i|j . 
The fact that the damping coefficient is twice the damping rate was explained in detail in [||. For the O(N) 
model with N = 4, 

2 rp2 



o(4) = 9A r 
' 16 



\J-Li 2 (e-™/ T ) , (56) 

7T TO V / 



in agreement with the result (4) of H. Note that in the large- iV limit, the damping coefficient vanishes 
- 1/N. 



V. THE 0(N) MODEL WITH SPONTANEOUSLY BROKEN SYMMETRY 

The O(N) symmetry of the Lagrangian j3^ ) is spontaneously broken to 0{N — 1) by taking to 2 < 0. 
Then, the potential ( p3| ) assumes the well-known "Mexican hat" shape, with the chiral circle |0 vac | = 
f n = (-m 2 N/4X)^ 2 as global minimum. Adding a small explicitly symmetry breaking term H <pi to the 
Lagrangian ( p2| ) "tilts the hat" in (/^-direction, such that the vacuum state is <^ vac = (f v ,0), where now 
/ w = {-rr^N/AX) 1 / 2 x 2 cos[a/3]/V3, cos a = (H A78A)/(-to 2 A/12A) 3/2 . The 0{N - 1) symmetry is 
restored taking H — > such that cos[a/3] — > y/3/2. 

Let us introduce new fields a = <p\ — f„ and 7r = (02, <^3, • • • , <?V), and corresponding masses: 

2 _ 2 i 12 A/ 2 2 2 4 A/ 2 
= TO + . m,, = to + — j^— • (57) 

Obviously, without explicit symmetry breaking (H — 0), the pions are true Goldstone bosons, to^ = 0. With 
these definitions the Lagrangian (p2|) becomes (constant terms are omitted): 



C(a, tt) = ~ (d M <7 0"a - to 2 a 2 ) + l - (d^ ■ - to 2 tt ■ tt) - l^L CT ( ( r 2 + tt • tt) - A( ff 2 + ff . ff) 



2 



(58) 



Taking N = 4, this is the Lagrangian of the linear a model. The parameters to 2 , A, and H of the original 
Lagrangian ( j32f ) are related to the physical meson masses m a = 600 MeV, to^ = 139 MeV, and the pion 
decay constant f„ = 93 MeV via to 2 = -(to 2 — 3to 2 )/2, A = N(m 2 — m 2 )/(8/ 2 ), H = m\]^. There are 
two types of interaction vertices in the Lagrangian (pq), a three-particle vertex proportional to Xf^ and a 
four-particle vertex proportional to A. Note that, since f n ~ A -1 / 2 , the three-particle vertex is formally of 
order X 1 / 2 . 

The er and 7r meson fields have the Fourier representation (cf. eq. (J^) 

a(t, k) = J d 3 x e _i k x <r(t, x) , 7r(i, k) = / d 3 x k x 7r(t, x) , (59) 
and are decomposed into soft and hard modes as follows: 

/d 3 k ( rl 3 k 

— e* kx a(t,k)e(fc c -|k|) , E(t,x)=y — e lk - x a(i,k)e(|k|-£: c ) , (60a) 

/d 3 k /" d 3 k 

^3e lkx 7r(t,k)6(fc c -|k|) , n(t,x) = y _- e ^- x 7r(t,k)e(|k|-fc c ) . (60b) 
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(The use of small letters is from now on reserved for the soft fields, the hard fields are denoted by capital 
letters.) 

The interaction between soft and hard fields is given by the action: 

Si[a,^,TZ,U} = S ( i XM [(7,^,tv,U] + S^ } [(7,E,7r,U} , (61) 

where 

S { j XM [a^,7v,U] =-i^A / ([3a 2 +7r-7r] £ + 2 [<r + £] tt • II + a [3 £ 2 + II • II] ) , (62) 

and 

sf°[o-,I!,7r,II] = —jy J ( 4 [a 2 + tt • tt] [a £ + tt • n] + 2 [3 cr 2 + tt • tt] £ 2 + 2 [cr 2 + tt • tt] n • n 

+ 4 [2 cr e + tt • n] tt • n + 4 [cr £ + tt • n] [£ 2 + n • n] ) . (63) 

Therefore, 

' 8Sl[(T, E, 7T, n] 

5a(x) 
4A/ W 



iV 
4 A 



( 6 ff(ar) (S(ar)) ffl # + 2 n(x) ■ (Tl(x)) 9 , w + 3 (£ 2 (a;)) ff , s + (n(x) • H(a;)) ff ,» 



- — ( [3cr 2 (x) + n{x) ■ *(x)] (E(aj))a >if + 2^(2;) tt(x) ■ (n(s)) 9ljf + 2 7r(a;) ■ (E(x)IL(x))a,n 

+ a(x) [3 (£ 2 (z)) ff . s + (n(as) • n(a;)) ff , # ] + (£ 3 (z)) ff , s + (£(x) n(as) • Il(x))*,* ) , (64) 

while 

^ 7 [a,S 7r,n] \ = _8A^ / _ _ + _ (na(a;)) _ + ^(3)) ) 

( 2Tt a {x) [a(x) (E(a;))a,7f + 7r(x) • (n(i)) ff , # ] + [cr 2 (a;) + tt(x) • tt(x)] (H^a:))*,* 
+ 7f a (ar) [(£ 2 (z)) ff , # + (II(a;) • U(x))^] + 2 [a(x) (£(z) U a {x))^ + n(x) ■ {U(x) U a {x))^} 
+ (£ 2 (:r) U a (x))^ + (U(x) ■ U(x) Ua(x))*,* ) • (65) 



A. Interactions between soft and hard fields to first order in A 

In this subsection, the expectation values (-)o-,7f on the right-hand side of eqs. (Q) and ( |65| ) will be 
evaluated in perturbation theory. Since A = N(m 2 r — to 2 )/(8/ 2 ) ~ 20 for N = 4 and realistic values 
of the parameters m a , m w , f v , this is certainly not a controlled approximation scheme. Therefore, the 
following results have to be viewed only as the first, but necessary, step to estimate the influence of an 
unobserved background of hard modes in the classical equations of motion for the O(N) model in the phase 
of spontaneously broken symmetry. 

To determine the interaction terms TZ a and TZ- Ka in the classical equations of motion ( |3lj ) to 
first order in A, due to the overall factors of A /„■ ~ A 1 / 2 and A in eqs. ( |64| ) and ( |65| ) the ex- 
pectation values in the terms proportional to A f n have to be computed only to order A f n , while 
those in the terms proportional to A have to be computed only to order 1. More explicitly, 
(E(x))CT,7f, (n o (a;))o-,7f, (£ 2 (a;))er, s , (II 2 {x)) s and (E(x) Jl a {x)) s have to be computed up to order Xf n , 
while (£ 3 (a;))^ s , (T, 2 (x)U a {x)) s ^ : (£(x) II 2 {x)) S fi, (U a (x)U b (x)) lf ^, and (II a (x) II 2 (a;)) ff ^ are required to 
order 1. For the latter, this of course means that the average ( • ) s ,jt can be replaced by ( • )o, and as a 
consequence, all expectation values vanish except for (H a (x)Ilb(x))o = cv^II 2 (x))q. For the non-vanishing 
expectation values one obtains: 
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/ <\ A n 


u > \ x 


(Un(x))* * 


.4AjV 


/ d 4 v 






= D^l(0)- 


1 N J 


( ^ 


(n^(x)) ff , ff 


= D ( *i(0)- 


1 N J 


1 d 4 y 



d 4 y ( (z - y)l 2 - \D& (x - y)] 2 ) 6 «7( tf ) 



d 4 y ( D^(x-y) '- D^(x-y) ' ) 2a(y) 



(66a) 
(66b) 
(66c) 
(66d) 



and 

(E(x) IIa(a;))g-,7f 



4AA 
N 



d A y D { ; ) {x-y)D^ , {x-y)~D { ° , {x-y)D^ , {x~y) 2n a (y) . (66e) 



Again, as discussed above in the symmetric case, all expressions are real-valued. For the left-hand side of 
the equation of motion for the classical a field one therefore obtains: 



£ a (x) = - 



n + m l(T) + ^(a 2 (x)+ir(x)-n(x)) 



a(x) 



N 



Za 2 {x) + tt(x) ■ tt(x) + 3 D^UO) + (N - 1) (0) 



v 7 i=l 



>(*) 



where 



d 4 2/ 



4' 1 (a - !/) - } (* ~ 6 <t(x) [3 cr 2 (y) + 7r(y) • n(y)] 



/ " r n 

TP(x)=i d 4 y 4^-j)-4 r) («-l/)^)-*(y)^), 



Tj 3 )( 3 



d 4 y 18 



D^(x-y) - D^(x-y) 



+ 2(N-1)[ D { ;\x-y) - ^(a-y) 



and m^(T) is given by: 



m^(T) = m 2 + 



4A 
~N 



3D^l(Q) + (N-l)D^l(0) 



(67) 

(68a) 
(68b) 

(68c) 
(69) 



The terms and the thermal mass correction are graphically displayed in Fig. ^. The conventions are as 
before in Fig. |^. A dashed line represents a pion and a solid line a sigma. Thin solid lines do not have a 
correspondence in the classical equation of motion, they correspond to a factor Aa(x) multiplying £ a (x) in 
the phase of the reduced density matrix (^7|). A filled dotted vertex corresponds to a four-particle ver tex A , 
a filled square vertex to a three-particle vertex A f n . Graphs (a,b) correspond to the two terms in eq. (|68a| ), 
graph (c) to (|68h|) , graphs (d,e) to the two terms in (38c) and (f,g) to [m^(T) — m%\cr(x). 
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(a) (b) (c) 



o 



(d) 



(e) 




x 
(f) 



i i 

x 

(g) 



FIG. 3. (a,b) The two terms in A sum over a is implied, (c) The term T^ 2 \ A sum over a is implied. (d,e) 

The two terms in T„ . (f,g) Order A contributions to the thermal mass. 

A Fourier transformation with respect to x yields the classical equation of motion for the momentum 
mode functions <x(t,k) (|k| < k c ): 



a(t,k) + i^ 3D^(0) + (JV-l)£>^(0) ^^(k) 



N 



4A_/ T /■ d 3 p 

N 

4 A /" d 3 pd 3 q 



/ 6(fcc " |P|) 6(fcc " |k " Pl) [3 P) * ( *' k P) + P) ' k " P)] 

0(k c - |p|)0(fcc - |q|)6(fc c - |k - p - q|) p)a(t, q) + *(t, p) • * (t, q)] k p q) 



TV 7 (2tt) 6 
( i ^) 2 E^ ) (*>k) = ^(t,k). 



(70) 



Here, = [k 2 + mj^T)] 1 / 2 . An explicit calculation of the terms T„ % \t, k) in the linear harmonic 

approximation is referred to Appendix C. 
For the classical pion fields one obtains: 



where 



□ + m 2 (T) + (a 2 (x) + * (x) • tt(x) + 2 / w a(x)) 



-a{ X )+{^f)j:r^{ X ), (71) 







'd { °\x 


-y)-D { °\x-y) 


2n a (x) [3a 2 (y) + 


7r(y) 


My)] ! 




L d *» 


'd { ;\ x 


-y)-D < Z\x-y) 


4a(x)a(y)n a (y) , 








L ^ 


D^\x 


-y)D^(x-y)- 


D^\x-y) D%\x 


-v) 


4 7Ta(y) , 



and 



4A 



mi(T)=ml + j^ D#(0) + (JV +!)£>£> (0) 



(72a) 
(72b) 
(72c) 

(73) 
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At first glance, m w (T) seems to be the thermal mass of soft (classical) pionic excitations. This, however, 
cannot be true, since in the chiral limit m n — > 0, the pions are massless Goldstone bosons, even at finite 
temperature |21|, whereas m^(T) as given by eq. j7^ ) is finite (cf. eq. (|43|)). This apparent violation of 
Goldstone's theorem will be resolved below. 

The terms T^J and the thermal mass correction are graphically displayed in Fig. ||. The notation is the 
same as in Fig. ^. Graphs (a,b) correspond to the two terms in eq. (72a), graph (c) to (721), graph (d) to 
(f72cl) and (e,f) to [m%(T) - ml}n a (x). 




, b 



V-V 



y s 



(a) 



(b) 



(c) 



X N 




— a 



(d) 



(e) 



(f) 



FIG. 4. (a,b) The two terms in T$ . A sum over b is implied, (c) The term T^J . (d) The term T^>. (e,f) Order 
A contributions to the thermal mass. 



(2) 



-(3) 



Fourier transforming eq. ([71]) with respect to x yields the classical equation of motion for the mode 
functions 7T (i, k): 



7T a (t,k) 



8XU f d 3 p 



N 



J e(fcc " |p|) e(kc ~ |k ~ p|) * (t ' p) * a( *' k p) 



4 A f d 3 nd 3 a 

~W J (2tt)6 (fcc ~ |p|)e(/Cc ~ 1^1)°^- - |k - P - q|) [ff(t, p)ff(i, q) + *(t, p) • * (t, q)] 7f„(i, k - p - q) 



4A/ 7r 
iV 



2 3 



53rW(t,k) = ^ a (t,k) 



(74) 



Here, = [k 2 + rn^T)] 1 / 2 . An explicit calculation of the terms Tj$(t, k) in the linear harmonic 

approximation is referred to Appendix D. 

For the variances of the noise terms one obtains with eq. (|26cj) up to order A: 



{l8 a(x)a(y) [d& (x - y) + D<?\x - y) 
+ 2 *(a0 • Tf (y) [l>« (x-y) + D« (i - y) 
(a - y)l 2 + [4' J (a: - y)} ") + (N - 1) ( [flW (a - y)f + [z^ } (a - y) 



4 A / ff 

iV 



{6a(a;)7f a (y) [d£> (x - y) + D<?> (x - y) 
+ 2 7f a (x) a(y) (x - y) + (x - y)] } 

2 {2 7ffc(») (*-») + (1 - 2/)" 
+ 2 <5 afc a(x) (7(2/) ( x - y) + flW (x - y) 



(75a) 



(75b) 
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2 6, 



ah 



D { ° ] (x - y) D ( *> (x-y) + D { °> {x 



y)D<$\x- 



II) 



(75c) 



This confirms the existence of correlations between the noise fields £ CT and ^ a as well as between ^ 7I - Q and 
£,ir b , a 7^ b, as mentioned at the end of section III. There also exists an obvious graphical representation for 
the variances. The difference to the graphs of Fig. || and |J is that one thick line on the vertex at space-time 
point y is replaced by a thin line, corresponding to the second factor Aa(y) or An^iy) m the phase in eq. 
(p7|). Also, due to £)>(x) + D <(x) = 2ReD™ + (x) one has to take the reaZ instead of the imaginary part of 
the respective diagrams. 



B. Spatially homogeneous solutions of the classical equations of motion 



In the following, let us focus on the time evolution of the zero-momentum mode functions a(t, 0) and 
7f Q (i, 0). For the sake of simplicity, let us also take the limit k c — » 0, i.e., only spatially homogeneous field 
configurations are considered to be classical, and let us assume t t — > — oo, tt — > +oo, to facilitate Fourier 

transformations. Then, the explicit form of the functions Ta% a (t, k) given in Appendices C, D shows that 

rj#„ (t, o) = rj 2 2 (t, o) = o, while 



ri 3 )(t,o) = -2P | 



27T TOo- — 0J 

9^ ro) K,0) + {N -l)M ( r\mq,0) 
2m rr 



<r(t,0) 



dt ff(t, 0) 



and 



Tj 3 > (t, 0) = -4 P / £ ^ , a (t , 0) _ 4 } , 0) 



2 



(76) 



(77) 



where {to, k) is the Fourier transform of (x) = D>) (x) (x) — (x) D^p (x), i,j = a or tt, cf. 

eqs. ( pi b| ) and ( |C2b| ). Using (2tt) 3 <5^ 3 ^(q) = VS^ q, where V is the total 3-volume of the system, as well as 
V Jd 3 p/(2?r) 3 = J2u for trie momentum integrals, and defining a{t) = a(t,0)/V, n(t) = H(t,0)/V, £ a {t) = 
£a-(t,0)/V, £-K a (t) = £-x a (ij 0)/V, one derives from fl70l ) the classical equation of motion for a(t): 



[d? + ml(T)} a(t) + 



4AA 

N 
4A 



3D { °U0) + {N-l)D { ;h0) 



4AA 
N 



+ "jy [^W + *(*) ' *(*)] + ^ 5 * = &(*) . 



(78) 



with 

ro*(T)=m2(T)+2 



4A/^ 
N 



do; 9 A^^ CT) (a;, 0) + (N - 1) {u, 0) 

27r mo- — uj 



4A 



7T 



iV 

2_ f4A/ T . 



3Di'2(0) + (JV-l)Df2(0)' 



iV 



d£ 



i 



+ (JV-1) 



d£ 



1 



(79) 



where use has been made of eqs. ( |69] ) and (C2b), and divergent terms have been removed by renormaliza- 
tion. The apparent singularity in the last integral poses no problem, since it is integrable. The dissipation 
coefficient is (cf. Appendix E): 



4 A U Y 2 9^^K, 0) + (N - 1) (m CT , 0) 



i( 7r7r )/ 



iV 



4A/ 7r V JV-1 



2m„ 



AT 



4 m 



8Trm 



Z coth^ 
4T 



(80) 
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The dissipation coefficient corresponds to the imaginary part of the diagram Fig. y (e), where the incoming 
(outgoing) a particle is on-shell and at rest. Therefore, the dissipation occurs physically due to the decay of 
the a into two it's The temperature dependence of r\ a is shown in Fig. [|. Note that even at T = 0, the 
dissipation coefficient does not vanish. This is physically plausible, since even then a a can always decay into 
two 7r's. This means, however, that dissipation (and associated fluctuations) persist even in the absence of 
a heat bath. In that case, the fluctuations have to interpreted as quantum rather than thermal fluctuations. 

For T = 0, and for the parameters of the linear sigma model, r\ a = 591.45 MeV, which is on the or- 
der m a , and thus quite large. In the chiral limit, — > and at T = 0, the dissipation coefficient 
rj a — » 3m^/(327r/^) = 745.26 MeV, which is even larger. r) a increases with T because of Bose-Einstein 
enhancement of the final two-pion state at finite T. 



2000 




m„ = 139 MeV 
m = 



200 - 
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FIG. 5. The temperature dependence of the dissipation coefficient n a in the case m n = 139 MeV (solid) and in the 
chiral limit m n = (dotted). 



For the variance of the noise field £ CT (t) the so-called "white- noise" approximation is employed (cf . Appendix 
G), which is consistent with the linear harmonic approximation (|A5| ) that made the equation of motion for 
the a field local in time and led to the term rj a dtd(t) in eq. (|78|). The variance of the noise field £, a {t) 
becomes 



S(t — t ') rju m<j coth 



2T 



(81) 



where (•){ denotes the average with respect to the Gaussian measure (J29J) . In the high-temperature limit, 
coth[77i CT /2T] — > 2T/to ct , and the variance coincides with what is known from the classical fluctuation- 
dissipation relation, which is employed in most treatments of the subject 0. However, due to m a ^> T 
for the range of temperatures of interest, this limit is not really applicable. Moreover, it would predict 
that the fluctuations vanish at T = 0, while the dissipation (p0[) persists. This certainly contradicts the 
fluctuation-dissipation theorem. The more general expression (|81| ) resolves this apparent contradiction, since 
coth[m CT /2T] — > 1 for T — > 0. Physically, the fluctuations at T = are quantum fluctuations originating 
from the decay of the a into two tt's. 

For the equation of motion of the pionic mode 7f a (t) = 7f a (t, 0)/V one obtains from eq. (|74|): 

[of + ml(T)] 7f„(t) + ^ a(t) 7f (t) + ^ [a 2 {t) + n(t) ■ n(t)] n a (t) + d t n a (t) = ^ a (t) , (82) 



where 



m^(T) = m^(T) + 4 



4A/ T 
N 



27T m-rr — UJ 



4A 
~N 
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and (cf. Appendix F) 




2 m w 



4mJ 1 — exp[— m^/T] 



4tt m% 



1 — exp[— m^./2m 7r T] exp[(m 2 — 2m 2 .)/2m 7r T] — 1 



(84) 



Damping of pions arises due to the processes 7r 7r — > cr and 7r cr — * 7r, where one 7r and the tr in the incoming 
channel come from the heat bath of hard modes. The temperature dependence of rj^ is shown in Fig. O. 
Note that r]^ is small as compared to in the temperature range of interest. The reason is the large a mass 
which strongly suppresses the phase space for the processes ir ir — > a and 7r a — > it (remember that one pion 
is at rest). Obviously, vanishes at T = 0, because then there is no background of hard 7r's or a's. It also 
vanishes in the chiral limit m v — > for all temperatures: truly massless Goldstone bosons are not damped 
(at least to first order in A and for k = 0). Correspondingly, the associated noise £ Va has to vanish as well. 
This can be explicitly seen from the expression for the variance of £„. a (t) in white-noise approximation (cf. 
Appendix G): 



Sab S(t — tf) Tfa inn COth 



2T 



(85) 



Finally, it should be mentioned that for the spatially homogeneous solutions under consideration, cross 
correlations between different noise terms vanish (cf. Appendix G). 
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FIG. 6. The temperature dependence of the dissipation coefficient r) n . 



C. Static, spatially homogeneous solutions of the classical equations of motion 



The equations of motion (|78h and Q82J J have the following interesting consequence. Let us consider a small 
perturbation (5a, Sir) of the vacuum ground state (a, 7r) vac = and let us consider the equations of motion 
only to lowest order in A/ ff ~ A 1 / 2 . All terms of order A as well as higher order terms in 6a and Sir will 
be neglected. To this lowest order in A, the variances of the noise terms vanish, cf. eqs. (|75|), thus also 
£<x = £,ir a = 0, and the equations of motion have the solution 
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5a = 



AT rrii 



Sir = 



(86) 



This static solution corresponds to a constant, temperature- dependent shift of the ground state <r vac 
— > (T vac = (5tj(T). This shift is identical to the well-known change of the vacuum ground state at finite 
temperatures [Q3j. To see this, let us consider the chiral limit — » 0, where (for AT = 4) A = m^./(2/^), 

and let us assume that T <C such that D^l)(0) may be neglected in comparison to D++(0) in eq. (|86|). 
Then one obtains with eq. (p3[): 



6- 



/ (1- — 



(87) 



in agreement with eq. (30) of J23J . 

Let us now resolve the apparent contradiction of Goldstone's theorem mentioned earlier. It has to be 
shown that in the chiral limit — > 0, the mass parameter for the static and homogeneous solution of eq. 
( p2| ) vanishes. In the present perturbative treatment, which is accurate to order A, one may employ the 
static solution ( [B6| ) for the a and ir fields in the equation of motion ( |S2|) in all terms which are proportional 
to at least one power of A /„• ~ A 1 / 2 . As shown above, the noise and fluctuation terms for the it field vanish 
in the chiral limit. Thus, to order A, one is left with the equation of motion: 



ml(T)7r a = , 



(88) 



where 



ml(T) 



lim 



4A 

~w 



D 



W(0) + (JV + l)uW ( o)_^( 
16 A [du>M<r\o>,ti) 



■■'>D { °h<S) 



(N- 



N 



2tt tUtt — ijJ 



(89) 



In the limit 



0, the last integral can be easily calculated (cf. eq. (p3T)) to yield 



KIT) . 



»K(0) + (AT + 1) D<;1 (0) - m (3 »S (0) + (W - 1) D« (0) 



+ 



16 A/ 2 

N rrvl 



(^(o)-^t(o)) 



(90) 



In the chiral limit, the coupling constant A — > m%./(2 f%), and all contributions to the thermal pion mass 
cancel, m 2 (T) — > 0, which completes the proof that Goldstone's theorem remains valid. In other words, 
the vanishing of rh^jT) implies that there exist non-trivial massless, static, homogeneous solutions to the 
equation of motion (Bq), which are, of course, nothing but the A^ — 1 Goldstone bosons. Note that the above 
cancellation of terms is equivalent to the arguments presented in |2l[ . 



VI. NUMERICAL SOLUTIONS 



In this section, numerical solutions of the classical equations of motion (78) and (|82j) are presented, to 
assess whether DCC's can form in the presence of dissipation and fluctuation. First, note that the sign of 
the quantity 

[mf(t, T)] 2 = rhl(T) + a(t) + ^ [a 2 (t) + *(*) • *(t)] (91) 

determines the time evolution of the pion fields in the absence of fluctuations or dissipation. For 
[m^ ff (£,T)] > 0, the pion fields simply perform oscillations with a constant amplitude, while for 

[m^ ff (t, T)] 2 < their amplitude grows exponentially. This exponential growth leads to large amplitude 
oscillations and, in turn, to a large number of pions in a given charge state (say 113 = no). The characteristic 
probability to find the 7r field aligned in a certain direction in isospin space leads to a probability for the 
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ratio R of neutral to all pions of P(R) ~ \j \f~R. This characteristic behavior was suggested as experimental 
signature for the formation of DCC's ||. Note that [m° ff (i, T)l can be negative only if the second term in 
(|9l|), i.e., the a field, is large and negative. An explicit calculation confirms that the mass parameter to^(T), 
eq. (|8^), is an increasing function of T. Therefore, exponential growth of the pion fields (and thus formation 
of DCC's) is most likely (and fastest) at T = Q. The following considerations will therefore be restricted to 
the case of vanishing temperature. In that case, however, there is no heat bath. Dissipation and fluctuation 
arise solely from the decay a — > tttt. 

The equations of motion (|78|), ( I84 ) are solved with a standard fourth-order Runge-Kutta method. The 
time step width was chosen to be At = 0.002 x 2-K/m a . The value of the fluctuating field £ CT is a Gaussian 
random number with variance rj a m a /(V At) p4| . It is chosen at the beginning of each time step and kept 
fixed during the Runge-Kutta step. 
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FIG. 7. The potential U(a) (in units of /*) as function of a (in units of /„■). The solid line is for m, = 139 MeV 
and the dotted line for the chiral limit m-n = 0. 

In Fig. ^ the potential 

U{a) = — a + — a + — a (92) 

is shown for N = 4, m a = 600 MeV, and A = (m| — mJ)/(2/|), where = 139 MeV, / OT = 93 MeV 
(solid line) . This function represents a cut through the potential energy surface of the linear sigma model 
at 7r = 0, cf. eq. (|58|). The absolute minimum, corresponding to the ground state, is at = 0. There is 
another local minimum at ~ — 1.87/ ff , and a local maximum at ~ — 1.13/^. In the chiral limit, 
m w — * 0, cr*- 2 ) — > — 2/ ff , while er^ 3 ) — * (cf. dotted line). Trajectories which are likely candidates for DCC 
formation obviously start at a < , and small non-zero values of tv (this is necessary because otherwise 
7r = remains a solution throughout the system's evolution). 

A representative candidate is shown in Fig. |8[ In part (a), the time evolution of a and 7r fields are shown 
for the initial conditions a — —1.14/^, 7fi = 0.002/^, ^2 — —0.001 f n . 7T3 = 0.001 f w . The derivatives of 
the fields are taken to be zero initially. Dissipation and fluctuation terms have also been set to zero for 
the solution shown in Fig. ||. One observes that the a field first "rolls" towards the minimum of the 
potential (p2|). Since this minimum is unstable in the direction of the 7r fields, the a field "rolls" on towards 
the absolute minimum = 0. During this process the tv fields grow. This growth is characterized by 
negative values of [mf{T)\ 2 or of sgn { [mf{T)] 2 } \mf(T)\, as shown in part (b) of Fig. | At this level, 
the evolution is conservative, and the system continues to oscillate around the ground state, since there is 
no way to dissipate the "potential energy" associated with the chosen initial state. This potential energy is 
simply converted into kinetic energy and vice versa in the course of the evolution. This also explains why 
the system periodically reaches the vicinity of the unstable minimum . This changes once dissipation 
and fluctuation is taken into account (it would also change if we solved the equations of motion in expanding 
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geometries The large-amplitude fluctuations of the pion fields make the observation of DCC formation 

experimentally possible. 

initial values (in units of fj: a=-1 .14, 7^=0.002, 7c 2 =-0.001 , 7c 3 =0.001 

- O) a a A A A 
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FIG. 8. (a) Classical trajectories without dissipation and fluctuations for a (solid), #1 (dotted), 7T2 (dashed), and 
7T3 (dash-dotted), in units of /„-. Initial values are a — —1.14/^, ffi = 0.002 /„-, 7T2 = — 0.001 f n , #3 = 0.001 /„■. (b) 
The function sgn {[«i^- ff ] 2 } \m^\ (in units of 
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initial values (in units of f„): <r=-1 .14, 7^=0.002, 7c 2 =-0.001 , tc 3 =0.001 
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FIG. 9. (a) Trajectories with dissipation and fluctuations for a (solid), 7fi (dotted), 7T2 (dashed), and 773 
(dash-dotted), in units of f n . Initial values are the same as in Fig. 8, the volume of the system is V = 1000 fm 3 . (b) 
The function sgn {[fC ff ] 2 } |wi,^- ff | (in units of f n ). 

Fig. M (a) shows a sample trajectory including dissipation and fluctuation terms in a comparatively large 
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volume V = 1000 fm 3 with the same initial conditions as in Fig. ||. Dissipation and fluctuations damp the 
oscillations of the fields, and they settle into the ground state. The time scale for this to happen is, for the 
sample trajectory shown here, about 10 fm. During this time, the pion fields oscillate strongly. Therefore, 
DCC formation would still be observable, if the system were to decouple after the first 10 fm. 

This conclusion holds, however, only for the one particular trajectory shown in Fig. ^, i.e., for one particular 
choice for the time evolution of the randomly fluctuating forces For a different random sequence, the time 
scale for damping could be larger or smaller. In a strict sense, one would have to average over an ensemble of 
time evolutions for £ CT for a given set of initial conditions for the fields [Q . The result of such an averaging 
is, however, predictable: on the average the fluctuations tend to cancel, and the ensemble averages of a and 
7r fields are constant in time. The assumed values depend, however, on the strength of the fluctuations, i.e., 
since £ CT ~ l/y/V, on the volume of the system. In large volumes, the fluctuations are not strong enough to 
"kick" the a field too far out of the ground state. On the other hand, in small volumes the fluctuations may 
have enough strength to let the a field reach the unstable minimum . Therefore, the ensemble average 
of the a field decreases towards — /„■ as the volume of the system decreases (on the average, the a field tends 
to be anywhere between the two minima and u^ 2 \ i.e., on the average, close to — The ensemble 
averages of the pion fields, however, are always zero (since the potential U(o~,tt) is symmetric with respect 
tO 7T — > — 7r). 

In Fig. [n] the time evolution of the fields is shown for a small volume V = 10fm 3 , for the same initial 
conditions (and the same random sequence for £ ff ) as in the previous two figures. As discussed above, in 
this case fluctuations are large enough to drive the system out of the ground state and induce disorientation 
of the 7r fields. One could speculate that the formation of DCC's is facilitated in smaller volumes, i.e., it 
should be more likely to observe them in collisions of lighter ions, or even pp-collisions. This could also 
provide an explanation for the CENTAURO events jll| , where heavy ions are not likely to play any role as 
collision partner. To confirm this, however, a more detailed investigation in an expanding geometry, with a 
realistic evolution for the temperature in the collision, and including modes with finite k (to study domain 
formation) is necessary fl6| . 

The problem with large-scale fluctuations in small volumes is, however, that one does not necessarily need 
to first restore chiral symmetry to observe them. In Fig. [H] a time evolution is shown in a volume V = 10 fm 3 
for the initial condition a = 0, 7fi = 0.002 f % , iT2 = —0.001 f % , tt^ = 0.001 /„., i.e., the initial values of the it 
fields are the same as before, but the a field is taken to be zero. The strong fluctuations £ CT drive the a field 
out of the ground state, and the small initial perturbations of the pion fields are strongly enhanced over a 
time scale of ~ 10 fm to produce large amplitude oscillations (and thus DCC formation). Note that the a 
field is always subject to these volume-dependent fluctuations, but if one starts with the true ground state, 
er = 7f = 0, as initial condition, the pion fields are not affected and remain zero throughout the evolution 
of the system. The conclusion would be that, in the presence of strong fluctuations, DCC formation is not 
necessarily a signal for restoration of chiral symmetry, a small perturbation of the ground state in the 7r 
direction seems to suffice. However, it is likely that this at first glance rather interesting phenomenon is 
physically identical with (and thus indistiguishable from) ordinary fluctuations in finite volumes. 

Note that there is a minimum volume V m in below which the dissipation coefficient rj a becomes zero. In 
order to have r\ a > 0, the a at rest has to be able to decay into two pions with finite, but opposite momenta. 
The lowest non-zero momentum state for a particle in box volume V = L 3 is k = (tt/L, 0, 0). Energy 

conservation in the decay process requires 2 = m a , or L mm = 2tt/ \Jm 2 a — 4raJ, or Vmin — 12.68 fm 3 , 
i.e., the case V = 10 fm 3 considered above is just on the order of the physically possible minimum volume. 
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initial values (in units of f n ): c>=-1 .14, 7^=0.002, 7c 2 =-0.001 , 7c 3 =0.001 
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FIG. 10. (a) Trajectories with dissipation and fluctuations for a (solid), ni (dotted), n2 (dashed), and 7r3 
(dash-dotted), in units of f w . Initial values are the same as in Fig. 8, the volume of the system is V = 10 fm 3 . 
(b) The function sgn {[mf ] 2 } \mf\ (in units of jv). 



initial values (in units of fj : o=0, 7^=0.002, 7t 2 =-0.001 , 7c 3 =0.001 
(a) 
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FIG. 11. (a) Trajectories with dissipation and fluctuations for a (solid), 7fi (dotted), TV2 (dashed), and 7T3 
(dash-dotted), in units of f n . Initial values are a = 0, 7fi = 0.002 f n , 7T2 = —0.001/^, 7T-3 = 0.001/^, the vol- 
ume of the system is V = 10fm 3 . (b) The function sgn {[m| ff ] 2 } |m| ff | (in units of f v ). 
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APPENDIX A: EVALUATION OF THE INTERACTION TERMS IN THE O(N) SYMMETRIC 

CASE 



In this Appendix, the evaluation of the interaction terms 7^ (i, k) in linear harmonic approximation is 
presented for the O(N) model in the symmetric case. Let us first define the functions 



.Mi (a:) =D>(x) -D<(a:) , 
M 2 (x) = Dl(x) - D\{x) , 
M 3 (x) = Dl(x) - D%{x) 

(called "memory kernels" in M), such that with eqs. (|39|): 



Mi(w,k) = e(|k|-fe c 



2tt 
2E^ 



[6{u; - Ek) - 5{lo + E k )] 



(Ala) 
(Alb) 
(Ale) 



(A2a) 



7H 2 (^,k) = 2tt 



d 3 p 



e(|p|-fe c )e(|k- P |-fe c 



( 27r )3 Mr, -W V," r, -W 4£pEk p 

{ ( [1 + ™(£p)] [1 + n(^k-p)] - n{E p )n{E^ p ) ) [ ^(w - E p - £ k _ p ) - 6{u> + E p + £ k _ p ) 

' [1 + n(E p )} n{E k _ p ) - n{E p ) [1 + n(E k _ p )] ) [ 5(u - E p + £ k _ p ) - % + £ p - £ k _ p ) ] } , (A2b) 



M 3 (u; 7 k) = 2tt | e( |p| _ fcc) e(|q| - k c ) 9(|k - p - q| - fc c ) gg ^ ^ 

x { ( [1 + n{E p )] [1 + n(£ q )] [1 + n(£ k _ p _ q )] - n{E p ) n{E^) n(£ k _ p _ q ) ) 
x 8{lo -E p -Eq- £k- P - q ) - 6{uj + E p + E cl + £ k _ p _ q ) 
+ 3 ([1 + n(E p )} [1 + n(£ q )] n(£ k _ p _ q ) - n(£ p ) n(£ q ) [1 + n(£ k _ p _ q )] ) 

x <5(w - £ p - £ q + £ k _ p _ q ) - £(^ + £ p + £ q - £ k _ p _ q ) } . (A2c) 

Note that all functions fulfill the symmetry relation A4i(uj, k) = — Mi{— (J, k). Then, eqs. ( |is| ) can be written 
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/H 3 k 
^3e' ;k - x ^ 3) (^o,k), 

where (for ^ — * — oo): 
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/•OO 

X i dr e~' iwr ^(xq - r, p) (p c (x - r, k - p) , 
Jo 

r a poo 
y^(x ,k) = J —M 3 (w,V)Q{k c -M)i J dre- MT ^(x -r,k) . 



(A4a) 

(A4b) 
(A4c) 
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Further evaluation is possible by employing the so-called linear harmonic approximation |8j for the time 
evolution of the classical fields: 



<p(t - t, k) ~ cos(£ k T) $5(t, k) - Sm( ^ kT) d t <p(t, k) 
— — ii k — 



(A5) 



This eliminates the r integrals over the history of the classical fields, rendering the final expressions local in 
time. Using the symmetry of Mi under to — > — to and the relation 



dr e 



iP- 



1 



X — UJ 



(A6) 



where P denotes the principal value, one obtains: 

y^&k)^ J ^!Se(fc c -|p|)e(fc c -|q|)e(fc c -|k- P -q|; 



did 

2^ 



(27r) e 
Mi(w,k) 



f 



f 



1 



4 i?2 £3 



1 



2 E\ E 3 



da; 
2^ 

(111.' 
2tT 



Mi(w,k) 



Mi(w,k) 



£■1 + E 2 + -E3 — a; i?i + E2 — E 3 — to Ei — Ei — E 3 — to 
x if b (t, k - p - q) <g(i, p) • q) 
1 2 f 



£1 -\- E/2 H- -£3 — uj £1 -+- £2 — -E3 ^~ ^ -Ei — -E2 — -E3 ^~ ^ 
X ip b (t,k - p - q) ft <g(i, p) • ft <g(t, q) 
1 1 f 



Ei -\- E2 H~ £3 — _£i -+- £2 — -£3 ^~ E\ — £2 ~\~ -£3 — uj 

f 



£*! — £*2 — -£3 — UJ 



ft ip b (t, k - p - q) <p(t, p) • ft <p(t, q) 



f 



AE 3 
1 

~8Ih 



[Mi(Ei +E 2 + E 3 ,k) - Mi{Ei +E 2 - E 3 ,k) +M 1 (E 1 - E 2 + E 3 , k) - Mi{Ei -E 2 - E 3 , k)] 

x tp b (t, k - p - q) Cp(t, p) • ft <p(t, q) 

[Mi(Ei +E 2 + E 3 ,k) + 2Mi(Ei + E 2 - E 3 ,k) + Mi{Ei -E%- E 3 ,k)] 

x ft <£ 6 (t, k p q) p) • q) 



+ 



f 



8 £/i £2 -£^ 



[Mi(^i + £2 + £ 3 ,k) -2M 1 (E 1 +E 2 - E 3 ,k) + Mi{Ei - E 2 - E 3 ,k)] 



d t tfb(t,k-p- q) d t tp(t, p) ■ d t fit, q) j , 



(A7) 



where E\ = -E k -p- q , £2 = Ep, £-3 = -E-q- For iV = 1, this expression reduces (up to a constant prefactor 
due to the difference between and .Mi) to eq. (50) of S. Furthermore, 



^ c 2) (i,k)^ |^e( fcc -|p|)e(fc c ~|k-p| 



do; 
2^ 



M 2 (w,k) 



1 



1 



1 



2 _£i £/2 

1 



dw 
2^ 



M 2 (w,k) 



E\ -\- E2 — uj E\ — £/2 ~~ uj 

1 1 



E\ -\- E2 — uj Ei — E2 — uj 



fb(t, p) Vc(t,k-p) 

ft p) ft <p c (t,k-p) 



4^1 
f 

4~B^ 



[A* 2 (#i + E 2 ,k) + M 2 (Ei - E 2 , k)] p 6 (t, p) ft p c (t, k p) 
[M 2 (£i + Ea,k) - M 2 (E! - E 2 ,k)} d t p t (t, p) <p c (t,k- p) 



(A8) 



where £a = -Ek-p, £2 = Ep- For AT = 1, this is equivalent to eq. (49) of (There is a factor of 2 missing 
in front of the last term of that equation.) Finally, 
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^ 3) (i,k)^-e(fc c -|k|; 



. du M 3 (uj : k) A^ 3 (^k,k) 

P / = ^ (t,k)H 7TF O t ip a (t,k) 



2tt £ k -u r v ' 2£ k 

which is equivalent to eq. (48) of ||. Then, the final expressions for the functions 7^(t,k) read: 
d 3 pd 3 q 

6 



(2tt) 6 
d 3 p 



■ e(fc c - |p|) e(fc c - | q |) [f (*» p) • £(*> q) + 2 p) <*)] 



xyfft.k-p-q) , 



rj 2 >(t,k)~ J A P ¥ e(fc c -| P |)^[(iv + 4)^(*,p)5 i , c + 4^(t,p)5 ac ] y«(t,k- P ) 



ri 3 )(i,k)^2(7v + 2)^ 3 )(t,k) . 



(A9) 



(AlOa) 
(AlOb) 

(AlOc) 



APPENDIX B: THE DISSIPATION COEFFICIENT IN THE O(N) SYMMETRIC CASE 

In order to derive the dissipation coefficient one has to compute the function 3(771, 0). In the limit 
k c — > and for uj = to > 0, k = 0, one first observes that the 6 functions corr espo nding to the decay of one 
particle into three and the reverse reaction (the second and third line in eq. (A2c)) have no support in this 
kinematic range. The only remaining contribution comes from the scat tering of the particle at rest with a 
particle from the heat bath of hard modes (the last two lines in eq. ( A2c)). Since the particles from the heat 
bath are supposed to be in thermal equilibrium, n(x) — [e x ^ T — 1] — 1 , the reaction rates observe the detailed 
balance criterion, or in other words, with the energy-conserving S functions one can rewrite the contribution 
from the reverse reaction in terms of a factor e~^l T times the original reaction rate. Then, 



d 3 pd 3 q 

W 



1 



M 3 (m,0) = ^ (1 - exp[-m/T]) J ^ 

x [1 + n(E p )] [1 + n(£ q )] n(£ p+q ) S (£ p+q + m - E p - £ q ) 



(Bl) 



(The last S function in the last line of eq. (A2c) has also no support for uj = m > 0.) This integral is 

i2 The angular 



most easily evaluated as follows. Let us first define p = |p| , q = |q| , E x = \fx 1 + 
integration (involving the angle between p and q) is substituted by an integration over -Ep+ q , with the 
Jacobian d-E p + q /dcos(p, q) = pq/E p+cl . This allows for a simple evaluation of this integral with the help 
of the 8 function: 



M 3 (m,0) = -^ (l-exp[-m/T]) f° ^ [1 + n(E p )\ f° ^ [1 + n(E q 
x n(E p + E q -m)0 (E p+q + to - E p - E q ) 6 (E p + E q - m - E p ^ q ) . 
The functions are equivalent to the constraints 

pq> (E p - to) (E q — m) > —pq . 



)] 



(B2) 



(B3) 



The right inequality is trivially fulfilled (E p > to), while the left is also true for all values of p and q, since 
E p < (p + to) 2 . Therefore, the functions do not impose additional constraints on the p and q integrations 
and can be simply omitted. Introducing the dimcnsionless variables a = m/T, x = E p /T, y = E q /T, 
abbreviating N[x) = (e x — and noting that 



(1 + N(y)) N(x + y-a) = N(x - a) — In J_J y 



(B4) 



one arrives at 



Mz{m,0) = ^L (1-e- 

32 7T J 



dx [1 + N(x)] N{x - a) In 



1 



(B5) 
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Introducing t = e a , substituting u — e x , and then z = (u — f)/(l — u), one obtains the final result 

3T 2 



M 3 (m,0) 



where 



32 7T 3 

- 1 dz 



Li 2 (e~ m / T ) 



Li 2 (t) = - / — ln(l - zt) 
Jo z 



(B6) 
(B7) 



is the dilogarithm, or Spence's integral. 



APPENDIX C: THE INTERACTION TERMS IN THE EQUATION OF MOTION FOR THE a 

FIELD 



In analogy to eqs. (Al) let us define 



Mf{x)=Df{x)-D^{x) , 



M% j) (x) = £>W (x) (x) - (x) (x) , i,j = a 



->(j) 



0), 



The Fourier transforms are rather similar to those in eqs. (A2) 



M { ;\w,k) = e{\k\-k c 



2tt 



IE 



d(cj-E^)-S(cj + E^) 



M u) (^k) = 2^ / ^ e(|p|- Ac )e(|k-p|-A c )— J- 



p ^k-p 



x {([l+n^W) 



1 + n (£ 



"k-p 



{^)n(E^ p )) 



sIu-eU-e^-sL + eU + e^ 



1 + n (e^ 



(4*)) [l + n«)]) 



(«-4°+^pH(«+4 -^)]} 



Then, the terms 7^ (a;) in the classical equation of motion for the a field read: 

d 3 k 



TP(x) = 6a(x) J -^e ik -S^(x ,k) , 
T^(x)^4^a(x) J ^e^S^Hxo^) , 



where (for ti 
5«(x 0) k)E 



— oo): 

^M^k) J e(k c -\p\)e(k c -\k-p\ 



2vr 



X! I dr e iWT [3 cr(a;o — r, k — p) cr(a;o — t, p) + 7r(xo — r, k — p) ■ 7t(:eo — r, p)] 
Jo 

5( 2 )(x ,k)^ | ^M^k) | ^e^-ipije^-ik-pi) 

/"OO 

x i / dr e -4 " 7- 7f a (x — r, p) a(x — r, k — p) , 
Jo 



S {3 \x 0l k) = 2 



du> 
2^ 



9A^^ ) (w,k) + (iV-l)X^ ) (w,k) e(A c -|k|)i/ dre-^ T a(x -r,k) 

i Jo 



(Cla) 
(Clb) 

(C2a) 



(C2b) 

(C3a) 
(C3b) 

(C3c) 



(C4a) 

(C4b) 
(C4c) 
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In the linear harmonic approximation (A5) these terms become: 

s«(*,k)~ J e(fc c -|p|)e(fe c -|k- P |) 



do; 
2^ 



2E[ g) E { ° ] 



p/^M ff) (-,k) 



1 



1 



E^+E^-o; E\°>-E 2 
1 



_J_p/"^^) (a;)k) 



er(f, p) <r(t, k - p) 

d t a(t, p) d t a(t, k - p) 



+ - u E^ - E^ - lo 



7r(t,p) • 7r(t,k- p) 



2E 



(°0 



X + , k ) - M? '[ E\°> - , k ) d t a(t,p)a(t,k-p) 



(<0 ( rp(cr) „(<r) 



d t Tr(t, p) • d t 7r(t, k - p) 



2E. 



(tt) 
2 

d 3 p 



^ K'+^^l-KM^-^k) t 7r(i,p)-7r(i,k-p) }> , (C5a) 



5( 2 )(t,k)^ | AP e(fc c -|p|)6(fe c -|k-p|) 
x{-ip/^W) 



+ - uj Ef> - E. 



(<0 _ pW 



2E^£ 2 



1 pfjl^k) 



TT a {t,p)a(t,k- p) 

9f7f (*,p)9 t CT(t,k-p) 



4E 



(a 



4E. 



' r M^ } f^+^.k) +MP (e^ -E«,k)] 7f (t,p)fta(t,k-p) 
^ (^ ) +^ ) ,k)-M l) (^-4 T) ,k)] ft7r (i,p)a(t,k-p) 



(C5b) 



5( 3 )(i,k)~-2 6(fc c -|k|) 



,, / 9M^ ) (u;,k) + (jV^l).M^ ) (^k) 
27r E ( ° ] - u) 

Ot <r\t, k) 



2E k 

where E± — E^_ and E% — Ep . The final expressions for T^*\t, k) are: 



(C5c) 



T«(t,k)^6 | (fc c -|p|)a(t,p)5W(t,k-p) , 

rj 2 ) (t, k) ~ 4 1 9 (fc c - |pI) £ p) (*, k - p) 



(27T) 

rj 3 )(^,k)^5( 3 )(t,k) 



(C6a) 
(C6b) 
(C6c) 



Since M^\<jJ,0) = 0, cf. eq. (|C2a|) , it follows from relations @ and @ that for fc c 0, 7; (1) (t,0) = 



Tj 2) (t,0) ee 0. 
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APPENDIX D: THE INTERACTION TERMS IN THE EQUATION OF MOTION FOR THE tt 

FIELD 



The interaction terms in the equation of motion for the tt held are: 



T^(x)^2ft a (x) 



d 3 k 

(2tt) 3 
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ikx P (1) (^,k) , 



d 3 k 

(2tt) 3 



(2tt) 3 



where pW(x ,k) = ^W^o.k), and ^(xo, k) = 5^ 2) (i , k), cf. eqs. (|C4a|) and (|C4bl) , while 

/J />OG 
(^k)0(fc c -|k|)i jf dre-^Tfa^o-T.k). 

In linear harmonic approximation we therefore obtain: 

7£>(t,k)~2 1 _^P_e(fc c -|p|)^(t,p)5W(t,k-p) , 

rW(t,k)~4 1 ^e(fc c -|p|)^,p)5( 2 )(i,k- P ) , 

7T a (*ik) 



(27T) 

rj 3 )(t,k)~-4e(fc c -|k 



da; M^^k) s , + ^ , A^(££ ) ,k) 



2.E, 
-(l)/ 



St7r (i,k) 

Again, since M^ J (u, 0) = 0, cf. eq. ( |A2a| ), it follows for k c -> that 7^(i, 0) = Tj a 2) (f., 0) = 0. 



2t i^-u, 



(Dla) 
(Dlb) 
(Die) 

(D2) 

(D3a) 
(D3b) 

(D3c) 



APPENDIX E: THE DISSIPATION COEFFICIENT FOR THE a FIELD 

The dissipation coefficient r\ a in the classical equation of motion for a(t), eq. (f7q) is defined by 



4 A U Y 2 9-M^ CT) (m CT , 0) + (JV - 1) (m CT , 0) 



AT 



2m„ 



(El) 



where ( w >k) is given by eq. (C2b). M^" 1 (m a , 0) describes the d ecay of a cr at rest into two it's (and 

the reverse reaction, corresponding to the second and third line of eq. ( |C2b ) ), an d the absorption of a cf by 
a a at rest, producing a a (corresponding to the fourth and fifth line of eq. ( C2b| )). 

However, & a at rest cannot decay into two ex's by energy conservation (similarly, two ex's cannot annihilate 
to form one), and also the absorption of a cx by a a at rest, producing a cr, is also impossible. Therefore, 
■M,' 70 ^ (m a , 0) = 0. Mathematically, this is immediately obvious inspecting the arguments of the 8 functions 
in eq. ( C2b| ) with ui = m a , k = 0. 

Similar arguments prevent the absorption of a n by a a at rest, producing a tt, but the decay of a heavy 
a at rest into two light 7r's (and the corresponding reverse reaction of tttt — > a) is always possible (provided 
m„ > 277^). Therefore (utilizing detailed balance, and taking the limit k c — > 0), 

d 3 k 



M { 2 nn) {m a ,0) = 2tt (1 - exp[-m CT /T]) 



(2tt) 3 



2E, 



1 + nlE, 



S[m„-2E, 



(tt) 



1_ / 1 _4m| 1 + exp[-m g /2T] 
8vr V m| 1 - exp[-m (J /2T] 



(E2) 



Note that this expression is non-zero even at vanishing temperature. The dissipation coefficient becomes: 

(E3) 



V<7 = 



4A/ 7r V N-l 



N 



4 m 



8tt m„ 



1 f coth-^ 

771' 4 J 
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APPENDIX F: THE DISSIPATION COEFFICIENT FOR THE tt FIELD 



The dissipation coefficient r/^ in the classical equation of motion for 7T a (t), eq. (183) , is defined by 



4A/ 7r 
N 



2 m,r 



(Fl) 



where .M^O^jk) is the Fourier transform of M^" ec l- (JC2bj) . M. y ^\ m ^^ 0) describes the deca y of a 7r 
at rest into a 7r and a a (and the reverse reaction, corresponding to the second and third line of eq. ( C2b| )), 
and the absorption of a tt or a by a tt at rest, producing a u or a 7r, respectively (corresponding to the fourth 
and fifth line of eq. (|C2b|)). 

The decay of a tt into cr and tt and the reverse process is kinematically forbidden. The remaining contri- 
bution reads with detailed balance (and taking the limit k c — > 0): 



(o-tt) 



r d 3 i 

(m 7r ,0)=27r(l-e X p[-m 7r /T]) y — 



d 3 k 1 + n ( E k 



2E 



to 



2E 



(tt) 



<5 (-E 1 , 



(tt) 



to 



1 

8ir ml 



•l 



4m^ 1 — exp[— m n /T] 



1 — exp[— m^ r /2mT T T] exp[(m^ — 2m^)/2m 7r T} — 1 



(F2) 



The dissipation coefficient therefore becomes 

,2 / 



4A/„ 



47r m 3 



'1 



i _ exp[— rrijr/T] 1 
1 — exp[— /2m n T] exp[(m^ — 2rn^)/2m. n T\ — 1 



(F3) 



This expression vanishes at T — and in the limit — > 0. 



APPENDIX G: THE VARIANCES OF THE NOISE FIELDS 



The variances of the noise fields are given by eqs. (|75|). In order to evaluate them, let us define in analogy 
to eqs. (CI) the "noise kernels" 



N^{x) = D { ;\x)+D^\x) , 

A/" 2 fe) (a:) = [x) D® (x) + D® (x) (x) , i,j = a or tt 



(i) (ij) 

Their Fourier transforms are quite similar to those of M. \ and M. 2 • 



Mi l \io,k) = e(\k\~k c ) 



2n 



2R 



(i) 



l + 2n 



7V 2 (ii) (a;,k) = 2tt 



d 3 p 



(2tt) 3 
{([l + n(4°) 



e(j P | - k c ) e(|k - p| - 



1 



P "^k-p 



1 + n [E, 



V k-P 



1 + n [R 



*(«-^+^ p )+«(a, + ^-JS^ p )]} . 



(Gla) 
(Gib) 



(G2a) 



(G2b) 



Using the fact that the hard modes are distributed according to the Bose-Einstein distribution function 
n(E) = {e E l T - I)" 1 , one can show that (cf. eq. (66) of g) 
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A/fV,k) = M[ l) {uj,k) coth 
A/" 2 (ij) (w,k) = M { 2 j) (cj,k) coth 



2T 



2T 



(G3a) 
(G3b) 



For the Fourier transforms of X a b(x,y) one therefore derives (£j — > 
l arT (k 0l k; g , q) = J' d 4 x d 4 y e l (k o x -^+ qo yo-^y) y) 



-co, i/ — > +00): 



^Pe(fc c -|p|)e(fc c -|k + q- P |) 

(2?r) 4 



x [ 18 cr(p , P) M a) ( k o ~ Po, k - p) cr(A:o + <?o - Po, k + q - p) 
+ 2^7f a (p ,p) A^ (7r) (fc -po,k- P) ^(fco + 9o -Po,k + q - p) 
+ (27r) 4 5(fc + go)^ (3) (k + q) [9AA 2 (CTCT) (fc ,k) + (7V-l)AA 2 ( ^ ) (fc ,k)] } , (G4a) 



X CT7ra (fc ,k;g ,q) = 



dpo d 3 p 

(2tt)4 



0(fc c - |p|)0(*c- |k + q-p|) 



Ga(p ,p)J\f^\k -p ,k- p) 7r Q (fc + g -p ,k + q - p) 
+ 2 7r a (Po, P) A^ (7r) Oq - po, k - p) a(k Q + q - p , k + q - p) 



(G4b) 



2-jr a 7r b (fco,k;(Zo,q) — 



dp d 3 p 

(2tt)4 



0(fc c - |p|)6(fc c - |k + q-p|) 



2 TTa (po, P) Ai (fc -Po,k- p) 7r 6 (fc + q -p ,k + q - p) 
+ 2S ab a(p ,p) Ni*'(ko -p ,k- p)a{k + q Q - p , k + q - p) 
+ 2 «U (2tt) 4 5(fc + g ) <5 (3) (k + q) (fc , k) } . 



(G4c) 



The focus of interest are the spatially homogeneous noise terms, k = q = 0. For this case, all integrals 
vanish in these expressions due to the fact that N[ 1 \oj, k) is proportional to 0(|k| — k c ), cf. eq. ( p2a ). This 
has the further consequence that all cross correlations between the noise fields vanish, I ab <~ S ab . Using 
(2Tr) 3 S^(k + q) = VS^l one arrives at: 



(k ,0;q ,0) = V 2n S(k + q a ) 
Tir a * b (ko,0;q ,0) = 5 a bV2ir5{k Q + q ) 



V N 



V N 



9AA 2 (CTCT) (fc ,0) + (N - 1)^2^(^0,0) 



2A/f" ) (fc ,0) . 



(G5a) 
(G5b) 



For the variance of the noise fields £ CT (i) and £n a (t) one therefore obtains: 



(&w&( t ')> {S /&ft0)£ ' (t ' 0) 



V 



V 



dfcp dq 
(2nf 



; i(fc t+go*') 



dfc o r ifa(t-t'i 1 /4A/„ 



2tt 



V 



2V 



9AA 2 (CTCT) (fc ,0) + (N- l)]4™\ho,G) 



(G6a) 
(G6b) 
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Here, (-)f is the average with respect to the Gaussian measure (29). Further evaluation is simplified by 



approximating Af^ a) (k a ,0) ~ Af^ a) {m a , 0), A/f^fco, 0) ~ 7V 2 (w) (m CT , 0), A/f^fco, 0) ~ ^\m v ,0), 
i.e., taking the energy fco to be the on-shell energy (i.e., since k = for both a and tt fields, the mass) of the 
respective particle. This approximation is consistent with the linear harmonic approximation which puts 
the energy in the M .2 functions on-shell. The consequence is that the fco~mtegral can be performed, yielding 
with eqs. flG3l) , ©, and @: 



- - i ^ ??<7 COth 



2T 



(&r<> (* ) (*' ) ) « - ^ 5 ( f _ f ' ) 77 TO t Vir Cot h 



2T 



(G7a) 
(G7b) 



The <5 function corresponds to white noise. Therefore, this approximation will be called "white-noise" 
approximation. 
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